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ABSTRACT 

We obtained a deep 150 ksec Chandra observation of the optically selected cluster of galaxies, RCS 231 8+0034, 
to investigate the gas mass fraction in this system. Combining our deep Chandra observation with an archival 
50 ksec observation, we derive gas mass fractions of fgas = 0.06 ± .02 and 0.10 ± .02 within r25oo and rjoo, 
respectively. The gas mass fraction in RCS 2318+0034 within r^m is typical of X-ray selected clusters. Further 
study shows that the large scale properties of RCS 2318+0034, including the relations between gas mass, X- 
ray luminosity and gas temperature are also consistent with the observed correlations of X-ray selected clusters. 
However, the gas mass fraction within r25oo is less than most X-ray selected clusters, as previously reported. The 
deep Chandra image of RCS 2318+0034 shows that this system is presently undergoing a major merger which 
may have an impact on the inferred gas mass fraction within r25oo. 

Subject headings: galaxies:clusters:general - intergalactic medium - X-rays:galaxies:clusters - 
cosmology : observations 



1. introduction 

The current matter in rich clusters of galaxies was accumu- 
lated from regions that would presently span approximately 
10 Mpc and should comprise a fair sample of the Universe. It 
has been known for some time that the dominant component of 
baryonic matter in rich clusters is the hot X-ray emitting gas 

■ (David et al. 1990; White & Fabian 1995; David et al. 1995; 
; Jones & Forman 1999; Ettori et al. 2003; Allen et al. 2004; 
'. VikhHnin et al. 2006; Allen et al. 2008; Sun et al. 2009; 

■ VikhHnin et al. 2009). White et al. (1993) showed that the 
] observed high gas mass fraction in clusters has significant cos- 

■ mological consequences regarding the total mass density of the 

■ universe. The segregation of baryons between stars and hot gas 
. depends on both the radius within a group or cluster and the 

total mass of the system. The ratio of gas mass to stellar mass 
\ increases from approximately unity in the cores of groups up to 
. more than than a factor of five at large radii in the richest and 
• hottest clusters (David et al. 1995, Gonzales et al. 2007). In 
' addition, since the gas is more extended than the dark matter, 
= the gas mass fraction increases with radius within groups and 
; clusters (David et al. 1995, White & Fabian 1995, Allen et al. 
, 2004, Vikhlinin et al. 2006). The greater extent of the gas, rel- 
ative to the stellar and dark matter, probably arises from some 
form of non-gravitational heating (e.g., mechanical AGN heat- 
ing and supernova driven galactic winds) which breaks the self- 
similarity between groups and clusters (Ponman et al. 1999, 
Arnaud & Evrard 1999, Lloyd-Davies et al. 2000, Bialek et al. 
2001, Neumann & Arnaud 2001, Finoguenov et al. 2002, Voit 
et al. 2005, Kravtsov et al. 2005, Dave et al. 2008, Bode et al. 
2009). 

Most studies of nearby clusters find that the gas mass frac- 
tion approaches a nearly constant value of fgas ~ 13% at large 
radii in the hottest and most massive clusters (David et al. 1995, 
White & Fabian 1995, Allen et al. 2004, Vikhlinin et al. 2006, 
Sun et al. 2009). This is slightly less than the WMAP value of 
17%. Based on these findings, there have been several cosmo- 
logical studies whose results are based on the assumption that 
the gas mass fraction in rich clusters is independent of redshift 
(Ettori et al. 2003, Allen et al. 2004, Ferramacho & Blanchard 



2007, Ettori et al. 2009, Juett et al. 2010). All of these studies 
have found results consistent with the concordant ACDM cos- 
mological model. However, a recent analysis of Chandra ob- 
servations of 13 optically selected, moderate redshift clusters 
culled from the Red-Sequence Cluster Survey (RCS; Gladders 
& Yee 2005) by Hicks et al. (2008), reports a significant scat- 
ter in the gas mass fraction of clusters with average gas mass 
fractions of fgas = 4.8% at r25oo and 6.8% at rsoo (where is 
the radius within which the average mass density of the clus- 
ter is A times the critical density of the universe at the redshift 
of the cluster). If this result is more representative of clusters 
compared to previous studies of X-ray selected clusters, then 
it would make precision cosmological studies based on the as- 
sumption of a universal gas mass fraction in clusters problem- 
atic. The cluster with the lowest gas mass fraction in the Hicks 
et al. (2008) sample is RCS 2318+0034 with reported values of 
/gai=4.0% and 1.5% at r25oo and rsoo, respectively. Due to the 
cosmological consequences of this result, we obtained an ad- 
ditional 150 ksec Chandra observation of RCS 2318+0034 in 
AO-1 1 to better constrain the gas mass fraction in this system. 

This paper is organized as follows. Section 2 contain a de- 
tailed discussion about the Chandra data analysis. Section 3 
presents evidence for a recent merger in RCS 2318+0034 and 
§4 examines the X-ray surface brightness distribution. All spec- 
tral analysis of the Chandra data is presented in §5. Section 6 
gives the gas mass, total gravitating mass and gas mass fraction 
at ^"2500 and rsoo and a comparison of our work with previous es- 
timates is given in §7. Our main results and their cosmological 
imphcations are summarized in §8. 

2. data analysis 

RCS 2318+0034 was observed with the Chandra ACIS-S 
detector on May 3, 2005 (ObsID 4938) and with the ACIS- 
I detector on Aug. 7, 2009 (ObsID 11718). Both observa- 
tions were taken in very faint (VF) telemetry mode. All data 
were re-processed with CIAO 4.3 and CALDB 4.4.6. We used 
acis^jrocess_events to apply the latest CTI and gain corrections 
and filtered the VF mode data to reduce the charged particle 
background. During the ACIS-S observation, the 12, 13, SI, S2 
and S3 chips were turned on and the cluster was positioned on 
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Fig. 1 . — left: hard band (2.0-6.0 kev) adaptively smoothed image, and right: .soft band (0.5-2.0 keV) adaptively smoothed image of RCS 2318-1-0034. The cross 
mai'ks the location of the peak surface brightness of the main cluster in the hard band image. 



the S3 chip. During the ACIS-I observation, the 10, II, 12, 13 
and S3 chips were turned on and the cluster was positioned on 
the 13 chip. Light curves were generated in the 2.5-7.0 keV and 
9.0-12.0 keV energy bands for the BI chip SI during the ACIS- 
S observation and for the BI chip S3 during the ACIS-I observa- 
tion. All time intervals with background rates in these energy 
bands exceeding 20% of the quiescent background rate were 
excised from further analysis. The resulting cleaned exposure 
times for the ACIS-S and ACIS-I observations are 50,155sec 
and I43,275sec. The appropriate ACIS blank field background 
files (Epoch D for the ACIS-S observation and Epoch E for the 
ACIS-I observation) were extracted from the CALDB and re- 
projected onto the sky to match the two RCS 231 8-1-0034 obser- 
vations. The same VP background screening was applied to the 
background data sets by only including events with "status=0". 
Finally, the exposure times for each chip in the background data 
sets were adjusted to produce the same 9-12 keV count rates as 
in the two RCS 2318H-0034 observations. 

Molecular contamination has been building up on the cold 
ACIS filters since launch. This only affects the ACIS QE be- 
low 2 keV. Since the ACIS-I observation of RCS 2318H-0034 
was taken more recently than the observations used to com- 
pile the background files, the throughput of ACIS at low en- 
ergies in the RCS 2318-1-0034 observation is less than it was 
during the background observations. To correct for this affect, 
we extracted spectra from a source free region on 13 beyond the 
virial radius of the cluster in the ACIS-I observation of RCS 
2318H-0034 and the 13 blank field background file. We found 
that the count rate in the 0.3-6.0 keV energy band (i.e., the en- 
ergy band we use to generate surface brightness profiles) dur- 
ing the RCS 2318H-0034 observation is still 4.6% less than the 
count rate in the 13 blank field data set even after adjusting the 
exposure time in the background data set to produce the same 
9-12 keV count rate. We therefore reduced the count rate in the 
13 background data set by 4.6% before subtracting the back- 
ground. We repeated this procedure for all 5 chips (all 4 ACIS- 
I chips from the ACIS-I observation and the S3 chip from the 
ACIS-S observation) used to generate background-subtracted 
surface brightness profiles for RCS 2318-1-0034. For all spec- 
tral analysis, we used the CIAO task specextract which extracts 
source and background spectra and computes photon-weighted 



effective area and response files. 

Due to the lack of measured galaxy redshifts for RCS 
2318-1-0034, we adopt a redshift of z=0.78 (which is the same 
as that used in Hicks et al. 2008) derived from the energy cen- 
troid of the H-like and He-like Fe-Ka emission lines in the 
ACIS spectrum. Assuming the standard ACDM cosmology 
with Ho=70 km s^^ Mpc^', £1^=0.3 and fiA=0.7, gives a lu- 
minosity distance of Dl = 4,870 Mpc to RCS 2318H-0034 and 
1" = 7.45 kpc in the rest frame of the cluster 



3. X-RAY MORPHOLOGY 

Adaptively smoothed ACIS images of the co-added RCS 
2318H-0034 observations in a hard (2.0-6.0 keV) and soft (0.5- 
2.0 keV) energy band are shown in Figure 1. In the hard band 
image, RCS 2318-1-0034 appears to be a fairly relaxed hot clus- 
ter, while the soft band image reveals that RCS 2318-1-0034 is 
undergoing a merger with a smaller and cooler system that is 
presently located about 100 kpc to the west of the dominant 
cluster The "x" in both figures indicates the location of the 
peak in the hard band surface brightness of the main cluster. 
The soft energy band image shows that the cool gas in the cen- 
ter of the main cluster has been displaced toward the east during 
the merger, possibly due to merger-induced sloshing (Ascasibar 
& Markevitch 2006). 

Since there are insufficient counts in the ACIS data to gen- 
erate a temperature map, we generated a hardness ratio image 
by adaptively smoothing the raw 0.5-2.0 keV and 2.0-6.0 keV 
images with the same kernel and then divided the two images. 
The resulting hardness ratio image is shown in Figure 2 along 
with surface brightness contours derived from the soft band im- 
age shown in Figure 1. The hardness ratio image shows that 
there is an elongated region of softer X-ray emission associated 
with the merging cooler system. The coolest emission from 
the main cluster is slightly displaced toward the east from the 
cluster core. The hardest X-ray emission arises from the region 
between the two merging cores. The spectral analysis presented 
below shows that this hard X-ray emission probably originates 
from a merger induced shock. 
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Fig. 2. — Hardness ratio image along with surface brightness contours de- 
rived from the soft band image shown in Fig. 1. 



4. SURFACE BRIGHTNESS PROFILE 

Due to the on-going merger in RCS 231 8H-0034, we extracted 
surface brightness profiles using both full annuli and annuli ex- 
cluding a 120° pie slice toward the west. All surface bright- 
ness profiles were computed using the location of the cross in 
Figure 1 as the origin. Data from the ACIS-S and AClS-1 ob- 
servations were co-added and the normalized blank sky back- 
ground images were used for background subtraction. To de- 
termine the appropriate weighting for the instrument and expo- 
sure maps, we extracted a spectrum from the central 500 kpc 
of RCS 2318H-0034 and fit the spectrum to a single tempera- 
ture apec model with the absorption fixed at the galactic value 
{Ngai = 4.04 X lO^" cm-2) and obtained kT=7.5 (6.4-8.6) keV, 
an abundance of Z = 0.12(< 0.33) (based on the abundance ta- 
ble of Grevesse & Sauval 1998) and an emission-measure of 
£ = (3.21 ±0.16) X lO^^cm^^. All error bars are given at the 
1 a confidence level. Instrument and exposure maps weighted 
by the thermal emission from a 7.5 keV plasma were then gen- 
erated for all 5 CCDs used to calculate the surface brightness 
profiles. Figure 3 shows the resulting background-subtracted 
and exposure-corrected 0.5-6.0 keV surface brightness profile 
excluding the emission from a 120° pie slice toward the west. 
We refer to this as our standard surface brightness profile. 

The X-ray surface brightness profile of clusters is commonly 
fit with a j3 model, given by: 



L(r) = Eo 1 




2\ -3/5 + 1/2 



(1) 



The results of fitting a j3 model to the background-subtracted 
and exposure-corrected surface brightness profiles using full 
annuli are shown in Table 1. This Table shows that as larger 
radii are included in the fitting process, larger values of j3 and 
core radii, r^, are obtained. The results of fitting the surface 




Fig. 3. — Background-subtracted and exposure-corrected 0.3-6.0 keV sur- 
face brightness profile of RCS 2318+0034 excluding the emission from a 120° 
pie slice toward the west along with our standard j3 model (solid line), also 
shown are the ACIS-I background (horizontal solid line) and r25oo and ''500 
(vertical dashed lines). 

brightness profiles that exclude the emission from a 120° pie 
slice are shown in Table 2. This table shows that smaller values 
of j3 and rc are obtained when the emission from the merging 
system toward the west is excluded, however, there is still a 
trend of increasing /3 and when larger regions are included 
in the fitting process. The main difficulty in fitting the surface 
brightness profile is that the X-ray surface brightness of RCS 
2318-1-0034 drops below the 13 background around 350 kpc, 
which is only about 2 core radii from the center of the cluster 
(see Figure 3). It is therefore difficult to obtain a tight constraint 
on the asymptotic slope of the X-ray surface brightness. 

To explore the systematic uncertainties in the derived param- 
eters of the j3 model, we generated a surface brightness profile 
excluding the emission from a 120° pie slice, using an exposure 
map weighted by the thermal emission from a 4 keV thermal 
plasma, since the gas temperature in clusters typically decreases 
by a factor of 2 at large radii (Vikhlinin et al. 2006). Table 2 
shows that the derived values for the /3 model are consistent 
with the values obtained using an exposure map weighted by 
the thermal emission from a 7.5 keV plasma. We also varied 
the background level by ±5% to determine the sensitivity of 
the derived parameters on the estimated background. This exer- 
cise shows that if the surface brightness profile is only fit within 
the central 350 kpc, then the results are fairly insensitive to the 
adopted background. However, if the emission from larger radii 
is included in the analysis, the results are quite sensitive to the 
exact background level used in the analysis. 

5. SPECTRAL ANALYSIS 

To further investigate the on-going merger in the center of 
RCS 2318-1-0034, we extracted spectra from the two cores and 
the region between the two cores. The results of fitting a sin- 
gle temperature apec model with the hydrogen column density 
fixed at the galactic value are shown in Table 3. While the un- 
certainties on the gas temperatures are fairly large, due to the 
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Table 1 

Surface Brightness Profile - Full Annuli 



Fitted Region 








(kpc) 




(kpc) 


(ct cm^^ arcsec^^ sec^') 


0-250 


0.64±.ll 


136 ±28 


(4.53 ±0.30) X 10-» 


0-350 


0.74 ±.09 


159 ±23 


(4.42 ±0.26) X 10"* 


0-450 


0.85 ±.09 


183±21 


(4.29 ±0.23) X 10"* 



Notes: Results of fitting a /3 model within different radii (Fitted Region) to the surface brightness profile extracted from within full 
annuli. The nominal background rates and instrument maps weighted by the emission from a 7.9 keV thermal plasma were used in 
the analysis. 



Table 2 

Surface Brightness Profile - Excluding a 120° Pie Slice 



Fitted Region 


1^ 


rc 


^0 




(kpc) 




(kpc) 


(ct cm~^ arcsec"^ sec 


0-250 


0.58 ±.09 


112±27 


(4.63 ±0.41) X 10" 


-« 


0-350'^ 


0.70 ±.09 


142 ±23 


(4.42 ±0.33) X 10- 


-8 


0-450 


0.75 ±.08 


152±21 


(4.36 ±0.23) X 10- 


-8 


0-250" 


0.58 ± .09 


112±26 


(5.39±0.47)xl0- 


-8 


0-350" 


0.70 ±.09 


141 ±23 


(5.07 ±0.38) X 10- 


-8 


0-450" 


0.75 ± .08 


152±21 


(4.99 ±0.36) X 10- 


-8 


0-250" 


0.59 ± .09 


114±27 


(4.60 ±0.40) X 10" 


-« 


0-350^ 


0.73 ±.10 


146 ±25 


(4.39 ±0.33) X 10" 


-8 


0-450^^ 


0.80 ±.09 


161 ±23 


(4.30±0.31) X 10- 


-8 


0-250^^ 


0.57 ± .09 


110±26 


(4.65 ±0.41) X 10" 


-8 


0-350^^ 


0.68 ± .08 


137±22 


(4.46 ±0.34) X 10- 


-8 


0-450^' 


0.71 ±.07 


144 ±19 


(4.41 ±0.32) X 10- 


-8 



Notes: Results of fitting a j3 model within different radii (Fitted Region) to surface brightness profiles extracted from annuli that 
exclude the emission from a 120° pie slice toward the west. The top 3 rows give the results using instrument maps weighted by the 
emission from a 7.5 keV thermal plasma. " Results using instrument maps weighted by the emission from a 3 keV thermal plasma. 
* Results obtained by increasing the background rates by 5%. Results obtained by decreasing the background rates by 5%. ^ Our 
standard J3 model. 



Table 3 

Spectral Analysis of Selected Regions 



Region 




kT 






(keV) 


Western Core 


3.24 


(2.70-4.19) 


Shock 


11.5 


(7.27-21.6) 


Eastern Core 


7.18 


(4.36-14.7) 



Notes: Best- fit temperatures and la uncertainties. 



limited photon statistics, we find that the western core is cooler 
than the eastern core, as inferred from the hardness ratio map. 
There is also a region of hotter gas between the two cores. This 
is consistent with a merger scenario in which the hotter gas was 
produced by a merger shock. The large uncertainties in the tem- 
perature prevent a more detailed analysis of the merger dynam- 
ics. 

A temperature profile was generated by extracting the emis- 
sion from within 3 concentric annuli with approximately 1,000 
net counts per spectrum. The spectra were fit to an absorbed 
single temperature apec model with the hydrogen column den- 



sity fixed at the galactic value. The resulting temperature profile 
shows some evidence for a radial decrease in the gas tempera- 
ture (see Figure 4). Also shown in Figure 4 is the broken power- 
law temperature profile used to estimate the total gravitating 
mass in RCS 2318+0034 and the resulting values for ^2500 and 
^500 (see §6). In addition, we extracted a spectrum from the in- 
nermost region excluding the emission within a 120° pie slice 
toward the west (which excises the emission from the merging 
subcluster), but the resulting gas temperature is consistent with 
the inner temperature shown in Figure. 4. Since the emission 
from RCS 2318+0034 drops below the X-ray 
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Fig. 4.— Temperature Profile of RCS 2318+0034 including the IfJ eiTor 
bars. The solid line shows the broken power-law temperature profile adopted 
from Vikhlinin et al. (2005) nomialized by T^p^c- The dashed lines show the 
broken broken power-law temperature profile normaUzed by the 1 a limits on 
Tspec- The vertical dashed lines show r25oo and rjoo- 

background at approximately 350 kpc, it is difficult to constrain 
the gas temperature at larger radii. 

6. GAS MASS, TOTAL MASS AND GAS MASS FRACTION 

Deprojecting the j3 model surface brightness profile, assum- 
ing a temperature-independent emission-measure, produces a 
gas density profile of: 




-3/3/2 



(2) 



where po is the central gas density. Thermal plasma models in 
XSPEC are normalized by the emission-measure, defined as: 



£ = 



10 



-14 



An[DA{l+z)] 



(3) 



where is the angular distance to the source, tig is the electron 
number density and «// is the hydrogen number density. The 
central gas density for a given j3 model is determined by insert- 
ing the density profile given by eq. (2) into eq. (3), integrating 
over a cylindrical volume with a radius of 500 kpc and equating 
the result with the observed emission-measure within this re- 
gion. Using our standard j3 model ()3 = 0.70 and = 142 kpc) 
and the observed emission measure within 500 kpc, we ob- 
tain ne{0) = (8.04 ±0.18) x 10"^^ cm^l The statistical un- 
certainty in the emission-measure produces an uncertainty of 
only 2.3% in the central gas density, since the density scales as 
the square root of the emission-measure. If we re-fit the spec- 
trum extracted from the central 500 kpc with the temperature 
fixed at 4 or 10 keV, we obtain best-fit emission-measures of 
(3.47 ±0.22) X 10-"* cm-5 or (3.17 ±0.15) x 10"'* cm^^ re- 
spectively. Thus, any variations in the gas temperature within 
this region would probably not generate more than an additional 
systematic uncertainty of 2.1% in the central gas density. 
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Fig. 5. — Ratio of the cumulative gas mass mass for several different 
p models to the cumulative gas mass derived from our standard j3 model 
(p = 0.70 and = 142 kpc). The lines correspond to: solid red line (/3 = 0.64 
and rr=136 kpc), dashed red line (/3 = 0.85 and r(.=183 kpc), sold blue line 
(j3 = 0.58 and (-,.=112 kpc) and dashed blue line (j3 = 0.75 and rM52 kpc). 
The vertical dashed lines show r25oo ^nd r^QQ. 

The dominant uncertainty in the gas mass in RCS 231 8-1-0034 
is not due to uncertainties in the central gas density, but is pri- 
marily due to the variations in the )3 model parameters derived 
from fitting the X-ray surface brightness over different regions 
(see the range in values of j3 and in Tables 1 and 2). Figure 
5 shows the cumulative gas mass derived from several differ- 
ent j3 models relative to the cumulative gas mass derived from 
our standard j3 model. The j3 models used for this comparison 
are given in the figure caption and are chosen to sample the full 
range of possible gas masses. Figure 5 shows that the range of 
acceptable j3 models produces an uncertainty in the gas masses 
at ?"25oo and r^oo of about 5% and 7-10%, respectively. We thus 
assign an average combined uncertainty in the gas mass at these 
radii of 10%. 

The total gravitating mass of a spherically symmetric system 
in hydrostatic equilibrium is given by: 



M,oti< r) 



kTr f dlnp dlriT 



^nipG \ dim dim 



(4) 



where k is Boltzman's constant, /i = 0.6 is the mean mass per 
particle, nip is the proton mass and G is the gravitational con- 
stant. The primary uncertainty in computing the gravitating 
mass in RCS 2318-1-0034 is the lack of a detailed knowledge 
of the gas temperature profile. We therefore compute the grav- 
itating mass assuming two temperature profiles: 1) an isother- 
mal profile and 2) the broken power-law temperature profile in 
Vikhlinin et al. (2005). We use our standard )3 model for the 
gas density distribution in both cases. 

Extracting a spectrum in an annulus from 70kpc to IMpc in 
RCS 2318-1-0034 and fitting the spectrum to an absorbed apec 
model yields a temperature of rvpec=6.32 (5.05-8.18) keV (la 
uncertainty). The resulting gas mass, total gravitating mass and 
gas mass fraction at r25oo and rsoo assuming the gas is isother- 
mal are shown in Table 4. The values for r25oo and rsoo 
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Table 4 

Results Assuming the Gas is Isothermal 



Tspec 


12500 


rsoo 


Mgas(r2500) 


Mgas(r500) 


Mtot (12500 ) 


Mtot(r50o) 


fgas(r2500) 


fgas(r500) 


(keV) 


(Mpc) 


(Mpc) 


(IQI^Mq) 


(IO^Mq) 


(IQi^Mo) 


(IO'^Mq) 






5.05 


0.308 


0.745 


0.796 


2.82 


9.96 


28.2 


0.080 


0.10 


6.32 


0.352 


0.837 


0.992 


3.25 


14.8 


39.9 


0.067 


0.081 


8.18 


0.408 


0.955 


1.25 


3.81 


23.1 


59.3 


0.054 


0.064 



Notes: The gas mass (Mgas), total gravitating mass (Mtot) and gas mass fraction (fgas) within r25oo and rjoo derived from our standard 
/3 model and the assumption that the gas is isothermal. Masses are shown for the best-fit emission-weighted temperature between 
70kpc and IMpc (Tspec) and the la limits on 7}pec- 



Table 5 

Results Using a Broken Power-Law Temperate Profile 



Tspec 
(keV) 


12500 

(Mpc) 


rsoo 
(Mpc) 


Mgas(r2500) 

(IQi^Ms) 


Mgas (rsoo) 

(lO^Mg) 


Mtot(r250o) 
(IOI^Mq) 


Mtot (rsoo) 
(IQi^Me) 


fgas(r2S00) 


fgas (rsoo) 


5.05 


0.320 


0.682 


0.853 


2.53 


11.2 


21.6 


0.076 


0.12 


6.32 


0.364 


0.765 


1.05 


2.93 


16.4 


30.4 


0.064 


0.096 


8.18 


0.420 


0.872 


1 .30 


3.43 


25.1 


45.1 


0.052 


0.076 



Notes: The gas mass (Mgas), total gravitating mass (Mtot) and gas mass fraction (fgas) within r2soo and rsoo derived from our standard 
P model along with the temperature profile in Vikhlinin et al. (2005). Masses are shown for the the best-fit emission-weighted 
temperature (lipec) and the 1(7 Umits on Tspec ' 



are calculated from the total gravitating mass profile. Since 
the resulting values for rsoo in Table 4 are slightly less than 
IMpc, we also extracted spectra from several annuli with in- 
ner radii of 70kpc and outer radii of 700, 800 and 900kpc. 
Due to the limited photon statistics in the Chandra observation, 
the best-fit temperatures are all consistent within the uncertain- 
ties. Table 4 shows that assuming an isothermal J3 model for 
RCS 23 18H-0034 produces gas mass fractions of 0.07 ± .02 and 
0.08 ± .02 at r2soo and rsoo, respectively. The error bars take 
into account the 10% uncertainty in the gas mass noted above 
and the uncertainty in Tspec • The gas mass fraction is essen- 
tially independent of radius at large radii since, for an isother- 
mal j3 = 0.7 model, Pgas ^ r"^ ' and p,„, ^ r^^. 

Based on the analysis of Chandra observations of a sample 
of 13 rich clusters of galaxies, Vikhlinin et al. (2005) found 
that the temperature profile of these clusters could be fit with a 
broken power-law model, given by: 

> ri.07, 0.035 < r/ri8o < 0.125, 

^l^n |l.22-1.2r/ri8o, 0.125 < r/ngo < 0.6. 

where (r) is the emission- weighted temperature excluding the 
central 70kpc (essentially Tspec)- For clusters, rsoo ^ 0.6ri8o. 
With this substitution, the temperature profile in eq. (5) is very 
similar to that derived for a sample of groups by Sun et al. 
(2009). 

Gas masses, total gravitating masses and gas mass fractions 
within r2soo and rsoo are shown in Table 5 using our standard 
j3 model and the broken power-law temperature profile in eq. 
(5) normalized by Tspec- Since eq. (5) depends on rsoo, we iter- 
ated in rsoo until a self-consistent solution was found. The self- 
consistent temperature profiles are shown in Figure 4. Table 
5 shows that assuming a broken power-law temperature profile 
for RCS 2318+0034 produces gas mass fractions of 0.06 ± .02 
and 0.10 ± .02 at r25oo and rsoo, respectively. The error bars 



take into account the 10% uncertainty in the gas mass and the 
uncertainty in Tspec- As opposed to the isothermal model, the 
gas mass fraction increases with radius due to the declining gas 

temperature. 

While RCS 2318+0034 is presently undergoing a merger, 
Nagai et al. (2007) have shown by performing mock Chan- 
dra observations of a sample of simulated clusters that X-ray 
derived gas masses in both relaxed and unrelaxed clusters are 
accurate to within 6%. Including this additional uncertainty in 
our derived gas mass would not a have a significant affect on 
our estimate of the gas mass fraction. Nagai et al. (2007) also 
found that X-ray derived gravitating masses are biased low by 
about 5-20% with no trend in cluster mass or redshift. This 
suggests that all X-ray derived gas mass fraction estimates are 
biased high by this amount. 

7. COMPARISON WITH PREVIOUS ANALYSIS 

Our estimated gas mass fractions using the combined 
200 ksec of Chandra data are approximately 50% greater at 
^■2500 and significantly greater at rsoo compared with the val- 
ues derived from the earlier 50 ksec observation (Hicks et al. 
2008). Comparing the derived parameters for RCS 2318+0034 
in Hicks et al. with our results, we find that both studies pro- 
duced similar values for the gas temperature (Tspec = 5.05 — 
8.18 keV compared to Tspec = 5-20 — 7.40 keV) and core ra- 
dius (rc = 142 ± 23 kpc compared to r^ = 171 ± 6 kpc) for the 
surface brightness profile. However, we obtain a value for the 
outer slope of the surface brightness profile corresponding to 
j3 = 0.70 ± .09 compared to j3 = 0.86 ± .03 in Hicks et al. This 
difference is primarily due to the poorer photon statistics in the 
initial 50 ksec Chandra observation. Larger values of j3 will 
over estimate the total gravitating mass and underestimate the 
gas mass and gas mass fraction. Jones & Forman (1999) found 
an average value of /3 = 0.6 when fitting the surface brightness 
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profile of a large sample of rich clusters observed with 
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gas mass fraction of fgas ('"soo) = 0. 10 ± 0.02. 
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Fig. 6. — Gas ma.ss fractions from Sun et al. (2009) (blue), Vikhlinin et al. 
(2006) (black), Allen et al. (2009) (green) and RCS 2318+0034 (red). Open 
symbols give fgu., within r25oo and solid symbols give fg^s within rsoo- 

the Einstein Observatory. More recently, in a study of 3 1 clus- 
ters culled from the REXCESS sample, which includes both 
relaxed and unrelaxed clusters, Croston et al. (2008) found an 
average power-law slope in the surface brightness profile be- 
tween 0.3r5oo and O.Srsoo among 6 keV clusters corresponding 
to J3 = 0.7. 

At r2500, we obtain Mgas{r25Qo) = 9.9l|^ x 10^^^^ ^nd 



M,o,(r25oo) = l-5lo 5 ^ ^^'"^ ^« assuming the gas is isother- 
mal. Hicks et al. (2008) derived MgaArisoo) = (9.7 ±0.6) x 
lO^^Mg andM,o,(r250o) = (2.43 ±0.28) x IQI^Mq assuming 
the gas is isothermal. However, the gas and total masses listed 
for RCS 2318H-0034 in Table 7 of Hicks et al. are inconsistent 
with the parameters given in Tables 3 and 4. The corrected val- 
ues are Mj,,„(r25oo) = 1 .60 x lO'^M© and M,„, (r250o) = 2.80 x 
lO'^ Mq (Hicks 2011, private communication). While there 
are some differences between our results due to differences in 
'"2500 and j3, our gas mass fraction at r25oo is consistent with the 
corrected Hicks (20 11, private communication) result. We find 
that assuming a broken power-law temperature profile increases 
both Mgas{r25(K)) ™d M,o,(r25oo) by 10% relative to the values 
derived by assuming the gas is isothermal. Thus, both temper- 
ature profiles produce the same value for the gas mass fraction 



of /ga.v('-2500) 

At rsoo. 



= 0.06 ±0.02. 

we obtain Mg„.,(r5oo) = 3.2; 



-1.0 
-0.7 
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M,„,(r5oo) = 4.0^;:^ X 10^^ Mq and ./;,„.,(r5oo = 0.06 ±0.02 
assuming the gas is isothermal. Hicks et al. (2008) de- 
rived Mg„,(r5oo) = (1.99 ±0.28) x IO^^Mq and M,„,(r5oo) = 
(12.9 ±2.0) X lO''* Mq. However, the corrected values are 
Mgasir50o) = 3.64 X IO^Mq and M,„,(r5oo) = 6.61 x lO^^ Mq 
(Hicks 2011, private communication). Assuming a broken 
power-law temperature profile reduces Mgasir^oo) by 10% (due 
to the smaller value of rsoo) and Mtotirsoo) by 25% (primar- 
ily due to the declining gas temperature). The assumption of 
a broken power-law temperature profile thus produces a larger 



Fig. 7. — Scaled gas mass vs. spectroscopic temperature for the X-ray se- 
lected clusters in Vikhlinin et al. (2006) (blue data points), the SZ-selected 
clusters in Andersson et al. (2010) (black data points) and the optically- 
selected cluster RCS 2318+0034 (red data point). 



8. DISCUSSION AND SUMMARY 

In general, the gas mass fraction in groups and clusters de- 
pends on the radius at which the gas mass fraction is computed 
and the virial mass of the cluster At small radii, the gas mass 
fraction also depends on whether the cluster has a cool-core. 
The central gas density in cool-core clusters is more than an 
order of magnitude greater than the central density in non-cool 
core systems. Since r25oo is only ^ 200 kpc in groups, the pres- 
ence of a cool-core can have a significant impact on the gas 
mass fraction at r25oo in groups. Sun et al. (2206) shows that 
there is indeed a strong anti-correlation between the gas mass 
fraction in groups at r25oo and the central gas entropy. For hot- 
ter systems (larger values of r25oo), the presence of a cool-core 
has less of an impact on the gas mass fraction at r25oo- 

Reports of the average gas mass fraction within r25oo vary 
from 4.5% in groups (Sun et al. 2009), up to 7.8% (VikhHnin 
et al. 2006) and 11.0% (Allen et al. 2008) in clusters. The re- 
ported gas mass fraction within r^oo varies from 8. 1 % in groups 
(Sun et al. 2009), up to 13.0% in clusters (VikhUnin et al. 2006, 
Ettori et al. 2010, Pratt et al. 2010, Juett et al. 2010). These re- 
sults are illustrated in Figure 6 and show that the differences are 
mainly due to the range in gas temperatures, or cluster masses, 
sampled in these different studies. The average gas temperature 
of the clusters in the Vikhlinin et al. sample is 4.5 keV compare 
to 7.7 keV in the sample of clusters analyzed by Allen et al., 
which contains a substantial fraction of clusters culled from the 
Massive Cluster Survey (MACS; EbeHng et al. 2001). 

Our estimate for the gas mass fraction within rjoo in the opti- 
cally selected cluster, RCS23 18-1-0034, is fully consistent with 
the results of X-ray selected clusters (see Figure 6). This re- 
sult is also illustrated in Figure 7, where we compare the scaled 
gas mass and spectroscopic temperature within rjoo for the X- 
ray selected clusters in Vikhlinin et al. (2006), the SZ-selected 
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clusters in Andersson et al. (2010) and the optically selected 
cluster, RCS 2318+0034. Croston et al. (2008) derived a 
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kT (keV) 

Fig. 8. — Bolometric X-ray luminosity vs. cluster temperature relation for 
the total emission within r^Qo. The RCS 2318+0034 data point is shown in 
red. All other data points are from Maughan et al. (2011). Relaxed clusters 
are shown in blue and unrelaxed clusters are shown in black. Also shown 
are the best fit relations for relaxed clusters (dashed blue line) and uni'elaxed 
clusters (dashed black line) derived in Maughan et al. 

scaling relation between the gas mass and temperature based 
on an analysis of clusters culled from the REXCESS sample. 
The scaling relation in Croston et al. is based on the temper- 
ature between O.lSrsoo and 0.75r5oo. Extracting a spectrum 
from within this region in RCS23 18+0034 gives a temperature 
of 5.3 ±091 keV. Using the scaling relations in Croston et al. 
gives an estimated gas mass at rgoo of (2.2 — 4.2) x lO'^M©, 
which is consistent with our estimate for the gas mass in 
RCS23 18+0034 (see Table 4). In addition, the scaling relation 
between fgas{r ^oo) ™d M,„f(r5oo) derived by Vikhlinin et al. 
(2009) predicts fgas{f 5qo)=^0.0%, which is in very good agree- 
ment with our result. Within r25oo, we find a gas mass fraction 
of 6% for RCS23 18+0034, which is significantly less than that 
found for most clusters (see Figure 6). Figures 1 and 2 show 



that RCS23 18+0034 is undergoing a merger which may affect 
the estimate of the gas mass fraction at r25oo. 

One of the most extensively studied cluster relations is that 
between the gas temperature and X-ray luminosity (e.g.. Edge 
& Stewart 1991, David et al. 1993; Fabian et al. 1994; Marke- 
vitch 1998, Arnaud & Evrard 1999). Early studies found a 
great deal of scatter in the Lx-T relation (Edge & Stewart 1991; 
David etal. 1993). Fabian etal. (1994) showed that most of this 
scatter was due to the inclusion of massive cooling flow clusters 
along with non cool-core clusters in the previous studies. The 
X-ray luminosity of massive cooling flow clusters can be up to 
five times greater than the X-ray luminosity of non cool-core 
clusters for a given temperature. By removing the emission 
from cool cores, Markevitch (1998) found that the scatter in the 
Lx-T relation can be greatly reduced. Similar conclusions have 
also been obtained by more recent studies (e.g., Maughan 2007, 
Prat et al. 2009, Maughan et al. 201 1). While most studies have 
shown that the Lx-T relation for clusters at a given redshift is 
not self-similar, there is significant evidence that the evolution 
of the Lx-T relation is nearly self-similar (e.g., Vikhlinin et al. 
2002, Maughan 2006, Maughan et al. 201 1). 

Maughan et al. (201 1) recently studied the Lx-T relation for 
a sample of 1 14 clusters observed by Chandra with redshifts 
of 0.1 < z < 1.3. The Lx-T relations in Maughan et al. were 
derived from both the total emission within rsoo and the emis- 
sion within (0.15 — 1.0)r5oo. Due to the limited photon statis- 
tics in the Chandra observation of RCS23 18+0034, the uncer- 
tainty in the gas temperature between (0.15 — 1.0)r5oo is too 
large for a useful comparison, so we only compare results us- 
ing the total emission within rsoo. Figure 8 shows the scaled 
Lx-T data for the relaxed and unrelaxed clusters in Maughan 
et al. (2011) along with the temperature and X-ray luminos- 
ity of RCS23 18+0034. Figure 8 shows that the temperature and 
X-ray luminosity of the z=0.78 cluster, RCS23 18+0034, is con- 
sistent with other unrelaxed clusters with similar redshifts given 
the observed scatter among clusters. 

In summary, we find that the large scale properties of the 
optically selected cluster RCS23 18+0034 are consistent with 
the large scale properties of X-ray selected clusters. At smaller 
radii, we find a lower gas mass fraction in RCS23 1 8+0034 com- 
pared to most X-ray selected clusters, which may be due to the 
on-going merger 

We would like to than A. Hicks for her help with compar- 
ing our results and also A. Vikhlinin for useful comments. This 
work was supported in part by the NASA grant G08-9122X. 
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